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Разработка алгоритма распознавания 
замкнутых плоских кривых” 


В данной работе продолжены начатые в [1-5] исследования по разработке алгоритма распознавания 
замкнутых плоских кривых (траекторий), построенных на двумерных сечениях Пуанкаре с помощью 
программы Модеег [1]. Предложенный алгоритм обработки изображений использует методику заливки и 
разработанный классификатор изображений. Алгоритм был реализован на языках высокого уровня 
С++и АсНоп ЗспрЕ 3.0 в среде ВоЙапа Ви|аег С++ 6.0 и Масготе а Е!азЬ АпипаНоп соответственно, 
протестирован на более чем 200 тестовых экземплярах изображений и показал достаточную эффективность. 


Введение 


Как известно [1], решение задачи прогнозирования регулярного и хаотического 
поведения нелинейных динамических систем сводится к анализу вида некоторых гра- 
фических изображений — сечений Пуанкаре. При этом сам процесс исследования можно 
разбить на несколько этапов. На первом этапе используется процесс дискретизации 
решения динамической системы с помощью методов численного интегрирования. На 
втором этапе проводится графическое и геометрическое моделирование полученных 
результатов. Если размерность системы не превышает трех, результат этого интегри- 
рования удобно представить в графическом виде на экране компьютера в виде после- 
довательности точек и образующих некоторую кривую (траекторию) в пространстве, 
характеризующую состояние системы в любой момент времени. На третьем этапе 
проводится анализ вида этой траектории и делается вывод о характере поведения сис- 
темы. В случае, когда размерность исследуемой системы превышает три, то результат 
численного моделирования системы в графическом виде в трехмерном пространстве 
можно представить только в виде сечений и проекций. Одним из видов таких сечений 
пространственных кривых является двухмерное сечение Пуанкаре. 


Постановка задачи 


Исходными данными являются изображения сечений Пуанкаре на сфере Пуас- 
сона, полученные с помощью разработанной программы Модеег [1]. Эти сечения 
построены в трехмерном пространстве и представляют собой сферу, на которой изоб- 
ражены множества точек (облака точек). В некоторых (регулярных) случаях эти мно- 
жества точек образовывают замкнутую трехмерную кривую — окружность, лежащую 
на поверхности сферы. Проецирование этой окружности на плоскость позволяет выде- 
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лить четыре случая: линия, окружность, эллиптическая и овальная кривая (рис. 1). Раз- 
работке алгоритмического и программного обеспечения для распознавания этих случаев 
посвящены дальнейшие исследования. 
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Рисунок 1 — Экземпляры распознаваемых изображений 


Таким образом, целью работы является разработка методов и алгоритмов реше- 
ния задачи распознавания замкнутых кривых, которые представляют собой траектории, 
построенные на двухмерных сечениях Пуанкаре. 

Объектом исследования выступают методы кластеризации и алгоритмы распозна- 
вания плоских кривых на двухмерных изображениях. 

Предметом исследования являются изображения специального вида с хаотичес- 
кими и нехаотическими траекториями. 


Алгоритм обработки и распознавания кривых 


В основе алгоритма лежит предварительная, основная и дополнительная обра- 
ботка изображения, которые проводятся в 6 этапов. Каждый этап получает на вход 


результат работы предыдущего шага и выдает результат для следующего. На рис. 2 
представлена схема взаимодействия модулей. 
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Рисунок 2 — Схема взаимодействия модулей 

Опишем алгоритм более подробно. 

Шаг 1. Предварительная обработка. Предварительная обработка необходима для 
устранения объектов на изображении, способных увеличить ошибку распознавания. 
Таковыми являются — вспомогательная сетка (фон), координатные оси, контуры объекта 
(опционально), одиночные пиксели, недостроенные траектории. 

Очистка изображения от вспомогательных линий и фона. Очистка производит- 
ся путем перебора всех пикселей входного изображения. При этом точки, имеющие 
черный цвет (#000000), заменяются на белый (# Е. Это необходимо для дальней- 
шего анализа траекторий, который будет использовать коды красного (#0000) и зеле- 
ного (#00100) цветов, описанного в шаге 4. 
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На изображении присутствует ось ординат (ОУ) красного цвета, являющаяся по- 
мехой для обнаружения кривых. Для глаза она имеет тот же оттенок, что и искомые 
кривые, но код её цвета не совпадает с ними. Ось располагается на заранее извест- 
ных координатах (199х0 — 199х299, 200х0 — 200х299, при разрешении изображения 
400х300 пикселей). Для корректного удаления оси нужно проверить соседние точки 
по горизонтали (координаты 198х0 — 198х299, 201х0 — 201х299) и произвести замену 
этими пикселями. Таким образом, удается избежать «рассечения» кривых. На рис. 3 
представлен алгоритм предварительной обработки изображения. 
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Рисунок 3 — Алгоритм предварительной обработки 


Шаг 2. Основная обработка изображения. В процессе построения кривых програм- 
мой-генератором возникают ситуации, при которых траектория является незавершен- 
ной. Эти случаи вызваны ограничением по времени, либо когда полное построение 
не представляется возможным. В виду того, что эта кривая строится по особому алго- 
ритму, представляется возможным восстановление потерянных (недостроенных) сег- 
ментов траектории. Образцы траекторий с различной степенью завершенности пред- 
ставлены на рис. 4 — 6. 
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Рисунок 4 — Траектории завершены на 50% 
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Рисунок 5 — Траектории завершены на 85% 
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Рисунок 6 — Траектории завершены на 100% 


Алгоритм восстановления отсутствующих сегментов траектории основан на цве- 
товом анализе изображения. Был применен пороговый метод - если яркость данного 
элемента превышает среднюю яркость локальной окрестности, тогда яркость данно- 
го элемента заменяется на среднюю яркость окрестности. Если полученная яркость 
превышает заданный порог, тогда точка принимает значение цвета кривой (#10000), 
в противном случае — белый (#1. 

Алгоритм эффективен в случаях, при которых потери сегментов составляют не 
более 20%. В остальных случаях появляется риск замены пустых областей внутри зам- 
кнутых траекторий на цвет контура (случаи, когда ширина замкнутой траектории в 2 
и более раз меньше длины). 

В зависимости от степени сложности (в это понятие входит процент потери сег- 
ментов на кривой, соотношение ширины к длине траектории, разброса и других по- 
мех) обрабатываемого изображения возможна повторная обработка. 

Шаг 3. Алгоритм распознавания кривых. Выделение областей изображения с из- 
вестными свойствами (признаками) кривых. 

На входе алгоритма имеется подготовленное изображение с объектами (#0000), 
очищенное от шумов, вспомогательной сетки и границ объекта. Так как исследуемые 
кривые были получены с помощью сечения трехмерной модели сферического тела, 
они обладают свойством симметричности по вертикали и располагаются строго сим- 
метрично относительно центра изображения. 

Исходя из этого, оптимальной точкой для начала поиска кривых является коор- 
дината 199х299 (для разрешения изображения 400х300). Направление поиска осуще- 
ствляется по формуле 199хМ (где М = 299...0) снизу вверх. Эта нижняя центральная 
позиция увеличивает вероятность обнаружения кривых, т.к. они обязательно пройдут 
через середину. Поиск осуществляется до тех пор, пока не выполняется одно из усло- 
вий — найден контур объекта, либо достигнута граница изображения. 

После того как был найден контур объекта, происходит обход вдоль контура по 
принципу «правой руки». Во время обхода контура оставляется «след» зелёного цве- 
та (#00100), нужный для остановки в случае движения по замкнутому контуру, а также 
визуального выделения области. Обход прекращается по достижении одного из усло- 
вий — следующий пиксель имеет зелёный цвет (#00100), либо достигнута граница 
сферы (на плоскости является окружностью). 

Контур считается выделенным. Продолжается поиск следующей траектории по 
тому же принципу, но уже начиная с координаты 199хК (где К — значение позиции пер- 
вого зеленого пикселя сверху по координате У). 
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Результат работы алгоритма, примененного относительно кривой, завершенной 
на 85%, приведен на рис. 7. 


Рисунок 7 — Траектории завершены алгоритмом 


Шаг 5. Классификация кривых. Поиск замкнутых, незамкнутых кривых и хаотич- 
ных областей алгоритмом заливки. 
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Рисунок 8 — Структурная схема работы алгоритма 


После выделения границ контуров происходит поиск траекторий, который зак- 
лючается в последовательном переборе всех пикселей изображения до нахождения 
пикселя красного цвета. Если пиксель был найден — создается временная копия изоб- 
ражения. На временной копии производится заливка найденной траектории синим цве- 
том (#0000ЕЕ), с последующим отсечением всех остальных траекторий по цветовому 
признаку. 
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После отсечения всех траекторий, кроме одной, производится определение ее раз- 
мера с последующей классификацией. Если линия слишком мала, то есть ее размер 
меньше порогового значения заданного, то она не подвергается дальнейшему анали- 
зу и считается шумом. В противном случае выбирается один из пикселей фона (цвет фо- 
на — #ЕЕЕЕЕР) и от него на временном изображении производится заливка зеленым 
цветом (#00ЕЕ00). После заливки производится поиск пикселей белого цвета, наличие 
которых свидетельствует о замкнутости траектории, а отсутствие о том, что она разом- 
кнута. Данный алгоритм адекватно работает лишь после подготовки изображения на пре- 
дыдущем этапе. На рис. 8 приведена структурная схема работы алгоритма. 

Шаг 6. На последнем этапе проводится дальнейший анализ изображения. Ана- 
лиз степени хаотичности распознанных кривых основывается на статистических дан- 
ных обработанного изображения. Анализируется количество замкнутых и незамкнутых 
траекторий, неправильно распознанных или несуществующих областей. 


Результаты тестовых испытаний 


Для тестовых испытаний были отобраны однотипные образцы траекторий. Цель 
испытаний — выявление случаев неэффективного либо неадекватного поиска и выде- 
ления кривых. Образцы представляли собой как завершенные, так и незавершенные 
объекты. Наибольшая эффективность алгоритма проявилась на изображениях с отсут- 
ствием хаотичных областей и завершенностью кривых на 85 — 100%. Ошибочное рас- 
познавание хаотичных областей проявлялось в случаях с плотностью расположения 
точек на расстоянии меньшем либо равным размеру самой точки. Результаты работы 
алгоритма распознавания приведены на рис. 9 — 11. 
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Рисунок 10 — Распознавание недостроенных траекторий 
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Рисунок 11 — Хаотическая область 


В ходе проведенного тестирования программа показала не плохую скорость 
работ, на картинках разрешением 400х300 пикселей анализ производился в течение 
1 — 3 с. Распознавание велось на картинках, смоделированных программой МОПРЕГЕК, 
при моделировании изображений были использованы параметры программы А = 1,70, 
В=0,9, С = 0,8, В = 1,91, & = 1, па1 = -0,8, пи2 = 0,6, пи3 = 0, шх = 1, шу = 0, 117 = 0. 
Выбор этих значений обусловлен тем, что при данных настройках возможно сгенери- 
ровать любой необходимый вариант траекторий для демонстрации возможностей ПО. 


Расчет сложности алгоритма 


Для оценки производительности алгоритма, его временной сложности, рацио- 
нально использовать асимптотический подход, т.к. алгоритм распознавания включает 
в себя более мелкие подгруппы алгоритмов и операций, зависящих от одних и тех же 
параметров и дающих несущественный прирост сложности при значительном увели- 
чении объема входных данных. 

За единицу времени принимается группа операций считывания одного пикселя. 
Характер алгоритма таков, что при увеличении количества траекторий на изображении 
количество анализируемых кластеров увеличивается, следовательно, увеличивается 
количество операций считывания/записи (заливка контуров). Характер сложности в 
данном случае носит следующую зависимость, которая представлена на рис. 12. 
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Рисунок 12 — График зависимости времени работы алгоритма 
от размера входных данных 
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Выводы 


В данной работе рассмотрены основные способы распознавания образов (изоб- 
ражений), а также их программная реализация. В ходе работы была произведена по- 
становка задачи распознавания замкнутых и не замкнутых траекторий на плоскости. 
В частности, был проведен анализ существующих методов распознавания изображе- 
ний, исследованы принципы и спецификации хаотических и нехаотических кривых 
на специального вида изображениях, были выделены основные признаки для построе- 
ния классификатора, разработано алгоритмическое обеспечение для распознавания 
кривых внутри класса, были проведены тестовые испытания и анализ результатов, 
произведен анализ эффективности алгоритма. 
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С.С. Литвин, К.А. Ручюн 

Розробка алгоритму розшзнавання замкнених плоских кривих 

У подашй робот! продовжен! розпочат! у [1-5] досллдження щодо розробки алгоритму розшзнавання 
замкнених плоских кривих (траекторй), побудованих на двовимрних перер1зах Пуанкаре при допомоз! 
програми Модеег [1]. Запропонований алгоритм обробки зображень використовуе методику заливки та 
розроблений класиф!катор зображень. Алгоритм був реалзований мовами високого р1вня С++ та АсНоп 
ЭспрЕЗ.0 у середовищ! ВоНапа ВиЙдег С++ 6.0 та Масготе а НазВ Апипайоп в1дпов1дно, протестований 
на б1льш н!ж 200 екземплярах для тестування зображень 1 показав достатню ефективн!сть. 


5.5. Гиут, К.А. Кис т 

Оеу@ортепЕ оЁ А!хогИ т о? Весосшйоп оЁ Фе СТозед Р1апе Сигуе5 

Те ргезепё \отК сопйпиез гезеагснез Бегип ш [1-5] оп 4еуеюртепе о# азогИ фт оЁ гесоот1 оп оЁ Фе 
с1озе р!апе сигуез (фа]есюме$) сопзасе4 оп Кмо-айпеп$!юпа| сиё оР Рошсаге Бу теапз оЁ ргоэтат 
Модеег [1]. ТЬе оЁеге4 а1еогИфит оЁ ипазе ргосеззше 115е$ а {есбаие оЁ роигпе ап4 Фе деуеореа 
апаПйег о{ ипазез. ТБе а]еогИфит Баз Бееп геа|зе4 ш 1апоиазез оЁ ШВ 1еуе| С +- ап4 Асйоп Зспр{ 3.0 ш 
фе епугоштепЕ ВоЙапа ВиЙ4ег С + 6.0 апа Масготефа Назв Апипайоп ассог@те]у, {ее оп тоге фап 
200 1е5{ сор1ез о{ пптазез ап Ваз зво\т за Йслеп еЁЙслепсу. 


Статья поступила в редакцию 25.10.2010. 
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